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Abstract 

We provide a tight approximate characterization of the n-dimensional product multicommodity flow (PMF) 
region for a wireless network of n nodes. Separate characterizations in terms of the spectral properties of appropriate 
network graphs are obtained in both an information theoretic sense and for a combinatorial interference model (e.g., 
Protocol model). These provide an inner approximation to the n 2 dimensional capacity region. These results answer 
the following questions which arise naturally from previous work: (a) What is the significance of l/v 7 "- m the 
scaling laws for the Protocol interference model obtained by Gupta and Kumar (2000)? (b) Can we obtain a tight 
approximation to the "maximum supportable flow" for node distributions more general than the geometric random 
distribution, traffic models other than randomly chosen source-destination pairs, and under very general assumptions 
on the channel fading model? 

We first establish that the random source-destination model is essentially a one-dimensional approximation to 
the capacity region, and a special case of product multi-commodity flow. For a wireline network (graph), a series of 
results starting from the result of Leighton and Rao (1988) relate the product multicommodity flow to the spectral 
(or cut) property of the graph. Building on these results, for a combinatorial interference model given by a network 
and a conflict graph, we relate the product multicommodity flow to the spectral properties of the underlying graphs 
resulting in computational upper and lower bounds. These results show that the 1/y/n scaling law obtained by 
Gupta and Kumar for a geometric random network can be explained in terms of the combinatorial properties of 
a geometric random network and the scaling law of the conductance of a grid graph. For the more interesting 
random fading model with additive white Gaussian noise (AWGN), we show that the scaling laws for PMF can 
again be tightly characterized by the spectral properties of appropriately defined graphs. As an implication, we 
obtain computationally efficient upper and lower bounds on the PMF for any wireless network with a guaranteed 
approximation factor. 

Index Terms 

Product multicommodity flow, wireless network, scaling law, capacity region. 

I. Introduction 

A. Prior Work 

An important open question in network information theory is that of characterizing the capacity region of a 
wireless network of n nodes, i.e., the set of all achievable rates between the n 2 pairs of nodes in terms of the joint 
statistics of the channels between these nodes. This has proved to be a very challenging question; even the capacity 
of a relay network comprising of three nodes is not known in complete generality. 

Instead of trying to characterize the capacity region for a general wireless network, the seminal paper by Gupta and 
Kumar [1] concentrated on obtaining the maximum achievable rate for a particular communication model, geometric 
random distribution of nodes, and randomly chosen source-destination pairs. They showed that the maximum rate 
for the protocol interference model scales as ®(\/^fn) for n nodes randomly placed on a sphere of unit area. This 
precise characterization has been followed by many interesting results for both combinatorial interference models 
and the random fading information theoretic model for large random networks; these include [2], [3], [4], [5], [6], 
[9] for communication theoretic models, and [7], [8] for information theoretic results. These results are crucially 
based on the assumption that a large number of nodes are randomly distributed in a certain region, and the inherent 
symmetry in the random source-destination pair traffic model. 

Since the relative locations of wireless nodes play an important role in the characterization of the capacity 
region, the notion of transport capacity was defined in [1]. A scaling law for the transport capacity for the protocol 
interference model was obtained in [1], while that for an information theoretic setting was obtained in [9]. The 
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transport capacity can be used to obtain an upper bound on the achievable rate-region for certain rate-tuples, but 
is not of much use in determining the feasibility of a certain rate-tuple. More recently, information theoretic outer 
bounds to the capacity region of a wireless network with a finite number of nodes were obtained in [10] for any 
wireless network, using the cut-set bound [11, Ch. 14]. We note that any achievable scheme can be used to obtain 
a set of lower bounds. While the above is only a discussion of a representative set of results in this area (see [12] 
for a more detailed summary), we note that there is no result which provides upper and lower bounds with a 
guaranteed approximation factor for a general wireless network with a generic random fading model. In this paper, 
we take the first steps towards providing such a tight characterization under very general assumptions. In doing 
so, we make connections between spectral graph theoretic results and network information theory. This results in 
efficient methods to compute the above tight upper and lower bounds. 

B. Contribution and Organization 

In Section [TTJ we consider the product multicommodity flow (PMF) as an n dimensional approximation of the 
n 2 dimensional capacity region. We show that the random souce-destination pair traffic model is a special case of 
PMF and it is essentially a one-dimensional approximation of the capacity region. 

In Section JIIJ we study the PMF for an arbitrary topology and a general combinatorial interference model, of 
which the protocol model is a special case. We show that the normalized cut capacity (equivalently conductance) 
of a capacitated network graph induced by the node placement and the interference model characterizes the PMF 
(within log n factor). For this model, we also obtain a precise scaling law for average delay using very elementary 
and almost structure independent arguments. We provide, possibly simpler, re-derivation of the (weaker by log 2 ' 5 n 
factor) lower bound on the maximum flow obtained by Gupta and Kumar for randomly chosen permutation flow on 
a geometric random graph with a protocol interference model. Our derivation illustrates the connections between 
the combinatorial properties of geometric random graphs and the maximum PMF. While we do not discuss in 
detail in this paper, the spectral properties of appropriate induced capacitated graphs characterize the scaling laws 
obtained for mobile networks in [2], [13]. 

In Section [TV] we address the question of characterizing the PMF for a wireless network with Gaussian channels 
and random fading. This is substantially more challenging than for the combinatorial interference model because 
there is no obvious underlying network graph that specifies the links which should be used for data transmission. 
We construct a capacitated graph whose cut capacity characterizes (in terms of tight upper and lower bounds) the 
PMF in the wireless network. This construction allows one to use classical network flow arguments to characterize 
and compute the PMF. We illustrate the generality of our results by obtaining scaling laws for a geometric random 
network and for a network where the number of commodities is constant. 

II. Traffic Flows 

In this section, we describe a class of traffic flows, namely, the class of product multicommodity flows, that we 
study in this paper, and its relevance. Consider a wireless network of n nodes and denote the node set as V = 
{1, ... , n}. A traffic matrix A = [Ay] G IR™ xn is said to be feasible, if for each pair of nodes 1 < i,j < n, 

data can be transmitted from node i to node j at rate A^. Note that whether a traffic matrix A is feasible or not 
depends on the model for the underlying wireless network, and we shall describe the precise models for wireless 
networks in the later sections. 

We denote the capacity region by A, i.e., A is the set of all feasible traffic matrices. Ideally, we would like to 
characterize A. However, this is a hard problem in most cases. Instead, we characterize an approximation of A 
under general assumptions on the wireless network. For this, we consider product multicommodity flow (PMF), 
defined as follows. 

Definition 1 (Product Multicommodity Flow (PMF)): Let node % be assigned a weight 7r(i), for 1 < i < n. 
Then the PMF corresponding to the weights it G MJf and a flow rate / G M + is given by the function [14] 

M : R n+1 i — ^ R nxn : 

tt(1>(2) ••• vr(l)vr(n) 

vr(2)7r(l) ... 7r(2)7r(n) 



M(/,7T) = / 



7r(n)7r(l) 7r(n)7r(2) 
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The PMF is an n-dimensional approximation of the n? dimensional capacity region A with product constraints. 
An important special case arises when all the weights are 1, i.e., = 1 for i = 1, . . . , n. We call such a flow 
uniform multicommodity flow (UMF). 

Definition 2 (Uniform Multicommodity Flow (UMF)): UMF with flow rate / G M+ is given by U(f) = fl, 
where 1 G M™ xn is a matrix with all entries equal to 1. 

We denote by /* as the supremum over the flow rates for which the PMF corresponding to the weights it is feasible, 
i.e., 

/; = sup{/Gi?+ : M(/,vr) is feasible}. 
We abuse notation and denote the corresponding quantity for UMF as simply /*. 

A. Inner Approximation to A 

We first show that the maximum UMF /* is a one-parameter approximation to the capacity region A. Consider 
the following parameter defined in terms of the capacity region as follows. 

Definition 3 (p*): For any A G M+ Xn , let p(X) = max, {££=1 X ik , Y2=i X ki} • Let 
L(x) = {A G K+ Xn : p{\) < x}. Then, define p* as 

p* = sup{x G M+ : L(x) C A}. 

Thus the quantity p* is a parametrization of a (regular) polyhedral inner approximation to the capacity region A. 
It is tight in the sense for any x > p* , there is an infeasible traffic matrix in the set L(x). 

Roughly speaking, the following result shows that UMF /* and p* are equally good approximations to the 
capacity region A. 

Lemma 1: If U(f) is feasible, then any A G R" xri such that p(X) < nf/2 is feasible. 

Proof: Consider any A such that p(A) < nf/2. Suppose that U(f) is feasible. Then there exists a transmission 
scheme such supports U(f). We now consider the two stage routing scheme of Valiant and Brenber [15] which 
routes U(p(\)/n) in each stage. Since U(f) is feasible, any A with p(X) < nf/2 is supportable by time sharing 
between the two transmission schemes corresponding to the two stages. To complete the proof of the Lemma, we 
need to describe this two stage routing scheme. 

In the first stage, each node i sends data to all the remaining nodes uniformly (ignoring its actual destination). 
Thus, node i sends data to any node j at rate Y^k x ik/ n < pW/ n - in the second stage, a node, say j, on receiving 
data (from the first stage) from any source i sends it to the appropriate destination. It is easy to see that due to 
the uniform spreading of data in the first stage, each node say j routes data at rate 2~2k x ki/n < p(X)/n to node 
i in the second stage. Thus, the traffic matrices routed in both the stages are dominated by U(p(X)/n). That is, 
the sum traffic matrix is dominated by U(2p(X)/n). Hence, if U(f) is feasible then p(X) < nf/2 is feasible. This 
completes the proof of Lemma [T] ■ 

Theorem 1: f* and p* are related as 

— < P 

Proof: Note that in general, the capacity region A may not be closed, and so we need a more careful argument 
We first show that < p*. By definition of sup it follows that for any e > 0, U(f* — ne/2) is feasible. Hence, 
from Lemma [T] any A G M™ x?1 such that p(A) < nf*/2 — e is feasible. Hence, again using the definition of sup, 
// "/• 2. 

Now for the other bound, assume that p* > nf*, and e = (p* — nf*)/2. Then, by definition of sup and p, 
U(nf* + e/2) is feasible, which is a contradiction. Hence, it follows that p* < nf*. 

m 

Thus, bounds on /* give bounds on p* which differ by at most a factor of 2. Subsequently, a scaling law for /* 
as a function of n is the same as a scaling law for p*, i.e. /* = Q(p*) as a function of n. 

'We present such a formal argument only once; similar arguments are implicit in many results that follow 
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The set of all feasible PMF clearly provides an n dimensional inner approximation to the capacity region, 
which is, in general, n 2 dimensional. Thus the characterization of the set of feasible PMFs provides a much better 
approximation to the capacity region than that the one-dimensional approximation given by set of feasible UMF. 
We next establish the equivalence of UMF with a traffic model with a randomly chosen permutation flow. 

1) UMF and Random Permutation Flow: In some previous work, (e.g., [1]), the capacity scaling laws were 
derived for the case where n distinct source-destination pairs are chosen at random such that each node is a source 
(destination) for exactly one destination (source) and such a pairing is done uniformly at random over all possible 
such pairings. Thus the traffic matrix corresponds to a randomly chosen permutation flow which is defined as 
follows. 



Definition 4 (Permutation Flow): Let S n denote the set of permutation matrices in I 
flow corresponding to a permutation E G S n and flow rate / G M + is given by S(f, E 



,nxn Then p ermutat j on 

) = /£. 



Many previous works study the scaling of /, where / is the supremum over the set of / G R+ such when a 
permuation E is randomly chosen from S n , the permutation flow S(f, E) is feasible with probability at least 
1 — 1/n 2 . We now show that when a permutation flow with flow rate nf and a randomly chosen permutation is 
feasible with a high enough probability, then the uniform mulicommodity flow U(f) can be "almost" supported 
when n is large enough. 



Lemma 2: For E G S n chosen uniformly at random, if (n/)E is feasible with probability at least 1 — ri 
0, then there exists a sequence of feasible rate matrices Y n such that 



-l-a 



, a > 



\Un(f) 



0(fn 



as n — » oo, 



where || • || denotes the standard 2-norm for matrice^l, and U n (f) is the uniform multicommodity flow for n nodes. 

Proof: From the hypothesis of the Lemma, it is clear that for at least (1— n~ l ~ a ) fraction of all n! permutations 
in S n , the permutation flow (n/)E is feasible. By definition and symmetry of permutations, we can write 

1 



U n (f) 



III 



Let us define the following indicator function 

h = 



1 (n/)Ej is supportable 
otherwise 

Consider a uniform time sharing scheme between all the re! permutation flows. Then the following traffic matrix is 
supportable. 

n! 

1 »i".nE, 



Thus 



\Un(f) 



re! ^-^ 



(a) 
< 



1=1 

n.\ 



j2\\(i-w(nm 



(b) 



X>-liKf < 



nf 



I n l+a 



n\ n 



0, as in oo. 



Step (a) uses triangle inequality for a norm and step (b) uses ||Sj|| = 1 for any permutation matrix Ej. 



Given a matrix M £ 



the 2-norm of M is ||M|| = sup{||Ma;|| : x £ R™, ||a;|| = 1}, where ||a;|| is the £2 norm of vector x G 
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From LemmaQ] if U (/) is feasible, then 5(n//2, X) is feasible for all S G 5 n . Thus, using an argument identical 
to that in the proof of Theorem [T] a scaling law for / is equivalent to a scaling law for /*, i.e. 

r = e(n/). 

Z?. Wireline Networks: PMF Over a Graph 

We briefly review the key results known for PMF on graphs with fixed edge capacities. These results will be 
useful in our analysis for PMF for wireless networks. 

Consider a directed graph G = (V,E), where an edge G E has a capacity C(i,j). Also, for ^ E, 
we take C(i,j) = 0. Then for a given it, f* for graph G = (V, E) is given by the solution of the following linear 
program (LP). 

maximize /, 

subject to ^2 (xijih k) - Xij(k,i)) = fir(i)iT(j), 1 < i,j < n, 

k:(i,k)eE 

y~] (xij(k,m) - Xij(m,k)\ = 0, VA; / l<i,j<n, 

m:(k,m)£E 
n n 

^J^Jscy(fc,m) < C(k,m), V(A;,m) G E, 
i=l j=l 

where the variables are / and {xij(k,m) : (k,m) G E,i,j, k,m = 1, ... , n}. The first two are flow conservation 
constraints and the third one is the capacity constraint. The total number of variables is less than 2n 4 and the total 
number of constraints is less than (n 3 + 2n 2 ). Hence, the above LP can be solved in polynomial time [16]. 

The well-known max-flow min-cut characterization for a single commodity flow naturally gives rise to the 
following question. Though the maximum PMF /* for a given weight vector can be computed in polynomial time, 
is there a corresponding result that relates /* and the properties of the graph. In their seminal paper, Leighton and 
Rao [17] obtained a characterization of /* in terms of the weighted min-cut of graph. We summarize their main 
result below. Let p n = \{i G V : ir(i) > 0}|, denote the number of nodes for which the corresponding element of 
7r is non-zero. Then, without loss of generality we assume that Yli=l 7r (*) = P*- 

Definition 5: For the graph G and weight vector n, define the min-cut by 

~ rr1 \ , J2(i,j):ieu,j&u c 

T(G, it) = mm — v ' , , — : , 

V ' ' ucv tt(U)tt(U c ) 

with notation that ir(S) = Y^ieS^i^ 1 ) ^ or an y set 

Theorem 2 (Theorem 17, [14]): In any directed graph G, the maximum PMF for weight it is related to T(G, it) 
as follows: 

V l0 gP7r / 

where the constants for the lower bound do not depend on the graph. 

Note that the upper bound follows easily because for a given PMF f w , the total flow from U to U c is 7r([/)vr(C/ c )/ 7r , 
which has to be less than the total capacity of the links from U to U c . The above characterization was crucial to 
the design of subsequent approximation algorithms for many NP-hard problems; a summary of these algorithms 
can be found in [14]. An important case of the above result is when ir(i) = 1 for all i = 1, . . . , n, i.e., the special 
case of uniform mulitcommodity flow. In this case, we have 

- T(G, 1) - mm , 

and 

n (Ji^Y\ < f* < t(GQ. 
Vlogny ~ J ~ v 1 
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III. Combinatorial Interference Model 

A combinatorial interference model defines constraints such that simultaneous data transmissions over only certain 
sets of links (or edges) can be successful. This is a simplified abstraction of a wireless network because in reality 
whether or not multiple simultaneous data transmissions are successful depends on the rate of data transmission 
and the interference power at the various receivers. We next describe the combinatorial interference model formally 
and illustrate it with example scenarios where this abstraction is a reasonable one. 

A. Model 

A combinatorial interference model for a given set of wireless nodes V = {1, . . . , n} defines the following two 
objects: 

(a) A directed graph G = (V, E) where E is the set of directed links (edges) over which data can be transmitted. 

(b) For each directed edge e G E, let 1(e) = {e£ E} be the set of edges (directed links) that interfere with a 
transmission on link e. Data can be successfully transmitted on link e at rate W(e) if and only if no transmission 
on any link in Z(e) takes place simultaneously. In general, the rate W(e) for a given power constraint can be 
different for different edges. The proof methods and results of the paper will not change (qualitatively) in this 
scenario. However, for ease of exposition we will assume W(e) = lo for all e£fi. 

We assume that for every edge (i,j) G E, edge G E, i.e., the graph G is essentially an undirected graph 
without the interference constraints given by the sets X(e)'s. This is a reasonable assumption in many time division 
and frequency division systems where the channels are reciprocal [18]. The interference sets l((i,j)) and T((j,i)) 
may not be identical because the transmissions which interfere with a signal received at node i may not be the 
same as transmissions which interfere with a signal received at node j. 

The above definitions can be used to induce a dual conflict graph as follows. 

Definition 6: The dual conflict graph is an undirected graph G D = (E, E D ) with vertex set E and edge set E D , 
where an edge e D G E D exists between e\ and e2 if e\ and e2 cannot transmit simultaneously due to interference 
constraints. Thus, each link e G E is connected to all links in 1(e). 

For the rest of the section, we will suppress the explicit dependence of all quantities on the combinatorial 
interference model parameterized by the graphs G and G D ; this helps to simplify notation. Let us denote the 
node degree and the chromatic numbei@ dual conflict graph G D by A = max eg £ |Z(e)| and k respectively. Note 
that k < (1 + A). Let {Ek},Ek C E, be the set of all possible link sets that can be active simultaneously, i.e., 
simultaneous transmissions on all the links in E^ at rate W(e) = 1 are feasible for the given interference model. 
Each Ek corresponds to a vector Ck G where Cfc(e) = lr eE E k } = l{e<=E k }- Let C be the convex-hull of all 
such vectors {C^}. Thus C is the set of all vectors C such that link capacities C(e) (for link e) can be obtained 
by time-sharing between the C^'s for the given interference model. We then define the capacity region to be the 
set of traffic matrices which can be routed over the graph G = (V, E), where each edge e has capacity C(e) for 
C G C. The formal definition is as follows. 

Definition 7 (Capacity Region (A)): The capacity region is the set of traffic matrices A G M nxn such that the 
following set of conditions are feasible for some C G C: 

^ (xij(i,k) - Xij(k,i)\ = Xij, l<i,j<n, 

k:(i,k)eE 

^2 (xij(k,m) - Xij(rn,k)j = 0, Vk^iJ, l<i,j < n, 

m:(fc,m)e_B 
n n 

^J^Jafy-(fc,m) < C(k,m), \/(k,m) G E, 
i=l j=l 

3 As long as W(e) is bounded below and above by a constant, scaling laws do not change even though the bounds for a given number of 
nodes n will change. 

4 The chromatic number of a graph is the minimum number of colors needed to color the vertices of the graph such that no two nodes of 
the graph which are connected by an edge share the same color. 
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where C(k,m) = C(e) with e = (k,m) G -E 1 ; variables are {xij(k,m) : (k,m) £ E,l < i,j,k,m < n}. 

Thus, the capacity region consists of all traffic matrices which are feasible using routing and link scheduling (time- 
sharing between the sets {E^}). We now illustrate this capacity region by a couple of special cases corresponding 
to widely used models for wireless networks. 

1) Protocol Model: The protocol model parameterized by the maximum radius of transmission, r, and the amount 
of acceptable interference, rj, is defined in [1] as follows. 

(a) A node i can transmit to any node j if the distance between i and j, rij is less than the transmission radius r. 

(b) For transmission from node i to j to be successful, no other node k within distance (1 + rj)rij (rj > a 
constant) of node j should transmit simultaneously. 

The corresponding definitions of E and Ed follow. A directed link from node i to node j is in E if < r. For 
a link, e G E, let e + denote the transmitter and let e~ denote the receiver. Then 

1(e) = {eeE: r s+e - < (1 + r/)r e+e -} . 

Thus the protocol model is a special case of the combinatorial interference model. 

2) SINR Threshold Model: Assume that all transmissions occur at power P, and the channel gain from the 
transmitter of node j to the receiver of node i is given by hij, i.e., if node j transmits at power P, the received 
signal power at node i will be hij P. A signal to interference and noise ratio (SINR) threshold model is parametrized 
by a threshold a such that a transmission from node % to node j is successful if and only if the SINR is above a, 
i.e., 

Pha 

3 > a 



T, k #Ph kj + N B 

For example, if we assume that each link transmits Gaussian signals and that the Shannon capacity on each link is 
achievable, then the threshold is given by N B(2 W - 1) (assuming W(e) = W for all e G E as before). 

We can define a corresponding combinatorial interference model such that the feasible simultaneous transmissions 
defined by the combinatorial interference model are a subset of that described by the SINR threshold model. Consider 
the set of directed links E 1 such that a link, from node % to node j, is in E 7 if and only if hji > 7. Also, 

define 1(e) = {e G E 1 : h e -e+ > 0}. Then link e can transmit at rate W if no other links in 1(e) transmit 
simultaneously, if and only if 7 and [3 are such that 

— - > a, Ve G £ 7 (2) 

It is easy to see that the above condition is satisfied if the following condition holds. 

B. Results 

We now derive results for the combinatorial interference model which relate the maximum PMF /* to spectral 
properties of the underlying graphs induced by the interference model. Most of the results in this section use ideas 
from known results. While important in their own right, these results and their proofs motivate the results for an 
information theoretic setting for wireless networks with Gaussian channels. Also, they provide alternate derivations 
for known capacity scaling laws in random networks. Towards the end of this section, we obtain simple bounds on 
the delay. 

1) Bounds on PMF: For any C G C, we denote the maximum PMF on graph G where each edge e has capacity 
(7(e), by /^(C) , and the corresponding min-cut by 

= mm — - ' J — — . 

v ' scv tt(S)tt(S c ) 

We denote the corresponding quantities for the special case of UMF by f(C) and *(C), respectively. Then we 
have the following lemmas. 
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Lemma 3: : C i— > R is a continuous function for any 7r > 0. 

Proof: Consider a cut 5 such that ir(S)ir(S c ) > 0. Then, the following is a continuous function of C: 

J2(i,j):i£S,j£S- C(.hj) 

tt{S)it(S c ) ' 

The lemma then follows since the minimum of a finite number of continuous functions is continuous. ■ 

Lemma 4: f n : C i— > E is a continuous function for any it > 0. 
Proof: For C G C and any e > define the set 

B S = {C £C: \\C-C\\ < 5}. 

To prove the lemma, we have to show that for any e > 0, there exists a S > such that for all C G 

|/x(C) : /,(C)|<e. 
For C G C consider 

<5i = min < aC(k, m) : C(k, m) > > , < a < 1, 

and 

5 = min < <5i , min 



i,j:7r(i)7r0)>0 2vr(i)7r(j) 

Then for any C G 65, it follows from[T]that f n (C) < f^iC) + e. It only remains to show that /^(C) > f w (C) — e. 
For this note that C >- C for all C G 65, where C_ is as follows: 



C(fc,m) 



C(k,m) = 

C{k,m) — 5\ otherwise 



Now by scaling all the variables by (1 — a) in the LP £[)) for C and using the monotonicity of /^(C) in C, we can 
see that f w (C) > (1 - a)f w (C) for all C G B 5 . If U{C) = 0, we are done. If not, choose a = min(e// 7r (C), 0.5), 
which gives /^(C) > f w (C) — e, and so we are done. 

■ 

We now define a quantity for the combinatorial interference model corresponding to the min-cut of a graph. 
Definition 8: The min-cut for the combinatorial interference model is defined as 

W _ = max mm — - — —— - — — — . 

n CaC ScV tt(S)tt(S c ) 

Note that ^* is well defined since ^(C) is a continuous function of C, and C is closed and bounded because 
it is the convex hull of a finite number of points. The above definition can be interpreted as the min-cut of the 
graph G, where each edge has capacity C(e), and the vector C is chosen from the set C such that it maximizes the 
min-cut of this graph G. The following result is an extension of Theorem [2] to combinatorial interference models. 

Theorem 3: f* is bounded as 

n (t^-) <K< %. 0) 

Proof: Since C is closed and bounded, it follows from Lemma[4]that there exists C* G C such that /* = f n (C*). 
Then, using Theorem |2j it follows that 

/* < tf(cr) < xi/;. 

Now, from Lemma [3l there is C G C such that = ty n (C). Using Theorem [2j it follows that 

f* > UC) = n — 1- 

VOgp n 

This completes the proof of Theorem [3] 
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Note that unlike the case for wireline networks (or equivalently graphs), /* is a hard quantity to compute. Also, 
note that ^ is a function of both G and the dual graph G D . We next relate the maximum UMF /* to spectral 
properties of graphs G and G D . 

Definition 9: The conductance of a graph G is defined as follows. 

*(G) = min E ^ Wcl[( *' j)g£1 , 
UCV,\U\<n/2 \U\ 

where In is the indicator function. 

Corollary 1: Recall that k is the chromatic number of the dual graph G D . Then, /* is related to &(G) as follows. 

n < r < 



Kn log n J ~ ~ n 

Proof: Consider vertex coloring for the dual graph G D = (E, E D ). The chromatic number of G D is defined 
to be k and hence we need k colors for vertex coloring of G D . Thus we have partitioned the set E into subsets, 
say, Ei, . . . , E K such that the links in each subset can transmit simultaneously at rate 1. Now let Cfc(e) = l{ eg £ fc }. 
Then, C corresponding to uniform time-sharing between the k edge sets E\ , . . . , Ek is given by 

C = -(C l + ... + C K ), 

K 

which is a convex combination of C\, . . . , C K G C. Hence, C(i,j) = 1/k for all i ^ j, and C € C. Then, using 
Theorem [2] and the definition of conductance above, 

/* > f*(c) > n (y—) = nf 



log n ) \ktl log n 

m 

For the upper bound, note that for any C G C, C ^ 1, i.e., C is lexicographically less than 1, and f(C\) > /(C2) 
if C\ > C2. Hence, /* = maxc e c /(C) < /(l). Then, the upper bound follows again by a straightforward use of 
Theorem |2] and the definition of conductance. 



2) Average Delay: We now provide bounds on the average delay for a class of traffic matrices. We measure 
delay in number of hops. We assume that the packet size is small enough so that the packet delay is essentially 
equal to the number of hops taken by the packet. This is similar to the assumptions in, for example, [1],[2], and 
[4]. 

We restrict ourselves to periodic link scheduling schemes (similar arguments extend to any ergodic scheduling 
scheme as well). For fixed networks, C is the convex hull of the set, {Ck}, which has a finite cardinality. Hence, 
any vector in C can be written as a linear combination of the C^'s. Thus to maximize the supportable uniform 
multicommodity flow it is sufficient to optimize over transmission schemes with periodic scheduling of links where 
the periodic schedule corresponds to time division between the C^'s. 

We obtain the following general scaling of delay. 

Theorem 4: Let S(n) be the total number of transmissions by the n wireless nodes on average per unit timeU 
When data is transmitted according to rate matrix A € A, the average delay, D(n), over all packets scales as 

D(n) = 6 , where A = £\. X tj . 

Proof: Let F denote the set of all possible paths (without cycles) in the network. The amount of flow generated 
at node i to be transmitted to node j is Ajj. Let us consider an arbitrary but nxec@ routing scheme where a fraction 
a]j of the flow from node i to node j is routed over path 7 G T. We assume that the traffic matrix A is feasible. 
Hence, there exists a link scheduling and routing scheme to support it. The total number of transmissions per unit 

5 The quantity S(n) is well defined since we consider periodic scheduling of links. 

6 Here, we consider a fixed deterministic scheme. However, it is easy to see that the result extends for any randomized scheme as well. 



10 



time at node I is Si j a lj^ij- Hence, the average number of transmissions per unit time in the entire network, 
denoted by S(n), is 

1=1 7 3i ij 7 i,j 

where if 7 is the number of hops on path 7. The total flow over a path 7 is J2i j a lj^ij> ie > tne fraction of total 
flow over path 7 is J2i j a lj^ij/~^- Hence, the average number of hops traversed by all packets is given by 



\^ ^ 13 J A 

7 1,3 



This completes the proof of Theorem |4] 



We note that the above result uses very little information about the specific underlying transmission scheme. 
Next, we present an immediate corollary of the above result and Theorem [J that bounds the delay for a scheme 
that maximizes the UMF in the network. 

Corollary 1: Since f* = ft f ^j^^ J. the corresponding delay scales as D(n) = O ( ^$^fy 

Proof: Consider the link scheduling scheme in the proof of Corollary [TJ where we partition the set of links E 
into subsets E\ , . . . , E K such that all the links in each subset E{ can transmit simultaneously. Note that this scheme 
can support UMF f = Q | ^^ n j, For this transmission scheme, every link transmits at rate 1 for at most 1/k 

fraction of the time. Hence, we have S(n) < Thus, it follows from Theorem |4] that D(n) = O ( ^|g£fr j ■ 



C. Computational Methods 

We now describe computational methods to obtain bounds on /* (the extensions to PMF are straightforward). As 
noted earlier, for wire-line networks, the computation of /* is equivalent to solving an LP. However, in a wireless 
network, the link capacity is a function of the link schedule. Since, the number of link schedules is combinatorial, it 
is a hard problem. Specifically, the question of checking feasibility of a rate vector A was proved to be NP-hard by 
Arikan [19], that is there exists an interference model and graph under which checking feasibility of A is NP-hard. 
Motivated by this, here we address the question of providing a simple computational method to bound /*. 

We use ideas of node coloring to induce a link schedule in a way similar to, for example, [20]. In particular, we 
can obtain an upper bound /* and a lower bound for maximum UMF /* in polynomial time such that 

/r < «/ 2 * 

The upper bound can be computed by solving the LP in (fl} with C(e) = 1 for all e£fi. For the lower bound, 
since the dual graph G D has chromatic number k, we can color the nodes of G D (which are given by the set E of 
wireless links) such no two nodes which share an edge share the same color. This in turn induces a link scheduling 
scheme, where each link in E is scheduled for at least a fraction 1/k of time, and the resutling C is such that 
C(e) > 1/k for all e £ E. Again, the lower bound can be computed by solving the LP in (fl} with C(e) = 1/k 
for all e G E. It is easy to see that /* < «/|. 

Now from Theorem |3l we know that Q, (^y^r^J < /* < ^- Thus, we can now also bound ^ as 

/*<*<0(/*logn) 

Thus, the upper and lower bounds differ by at most a factor of k log n. In addition, using the algorithm in [17], we can 

find a vector C(e) and the corresponding cut (U, U ) such that the capacity of this cut, min,s c y — ^^rfffpr > 

is within a factor K(logn) 2 of ^. 
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D. Application 

We now illustrate our results for the combinatorial interference model through an application to geometric random 
graphs. The geometric random graph has been widely used to model the topology of wireless networks after the 
work of Gupta and Kumar [1]. However, it has been a combinatorial object of interest for more than 60 years. We 
derive scaling laws for a combinatorial interference model which is more restrictive than the protocol model. Note 
that the lower bound hence obtained is also a lower bound for the protocol model. Specifically, the lower bound 
is weaker by log 2 ' 5 n compared to the lower bound obtained in [1]. Specifically, we show that the scaling of the 
lower bound is closely tied to the known combinatorial properties of geometric random graphs. 

We first define the restricted protocol interference model. It is also parameterized by the maximum radius of 
transmission, r, and the amount of acceptable interference, 77. 

(a) A node i can transmit to any node j if the distance between i and j, is less than the transmission radius r. 

(b) For transmission from node i to j to be successful, no other node k within distance (1 + 77) r (where 77 > is 
a constant) of node j should transmit simultaneously. 



We now state a version of the well-known Chernoff 's bound for binomial random variables that we use multiple 
times in this paper. 

Lemma 5: Let X\, . . . ,Xn be i.i.d. binary random variables with Pr(Xi = 1) = p. Let S n = X)fc=i f° r 
n = 1, . . . , N. Then, for any 5 G (0, 1) 

Pr (\S n — np\ > 5np) < 2 exp 



Specifically, for 5 = ^/ 2I ' n ° sn , we have 



Pr (\S n — np\ > yj 2Lnp log n^j < —. 



Consider n wireless nodes distributed uniformly at random in a unit square, and the interference model given 
by the protocol model with transmission radius r. We denote such a wireless network by G(n,r). It is well- 
known that for G(n,r) to be connected with high probability, it is necessary to have r = f2 (y 7 log n/n). We 
take r = 6 (log 3 / 4 n/y/n) and prove the following bounds on the maximum UMF, /*, for the restrictive protocol 
model; the lower bound is only log 2 5 n weaker than the result of Gupta and Kumar for the protocol model with 
r = B(y / logn/n). 

Lemma 6: For G(n,r), with r = 0(log 3//4 n/y/n), maximum UMF is bounded as 

n ( „ * 575 ) < r <o ( 1 



\ ti 3 / 2 log 5//2 n / V n 3 / 2 log 3//4 77 , 

Proof: To prove the above bounds, we obtain appropriate upper and lower bounds on the quantity ^. These 

bounds along with Theorem [3] imply Lemma [6] To obtain an upper bound on ^t, we evaluate the cut-capacity for 

a specific cut-set. For the lower bound, we first, establish that a grid graph on n nodes is a sub-graph of G(n, r) 

and then use the known conductance of the grid graph. 

First, consider the upper bound on <I>. Specifically, consider the square, say S, of area 1/9 (of side 1/3) that is 

in the center of the unit square. Let S be the set of nodes that fall inside this square. By definition, we have 

* < = sup f—^ . 

cec \o\\b c \ 

Therefore, it is sufficient to required obtain bound on ^f(S). 

Corresponding to node i, define a random variable Xi G {0, 1} which is 1 is i is in S, and otherwise. Since 
nodes are placed uniformly and independently at random in the unit area square, Xi are i.i.d. binary random variable 
with Pr(Xj = 1) = 1/9. Now, Yj£=i Xi is the number of nodes in S. Using Lemma [5] with 6 = 0.009, it follows 
that for large enough 77, \S\ G (O.ln, 0.2n) (and so |5 C | G [0.8n, 0.9n]) with probability at least 1 — n~ A . Now, 
consider squares S ^ 1 of sides 1/3 + 2r and 1/3 — 2r respectively with their centers being the same as that of 
S. That is, S 1 C S C S°. Let A = S° - S and A 1 = S - S\ Thus, .A 1 is a strip of width r surrounding S of 
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and .4° is a strip of width r on the boundary and inside S. Since, r = 0(log 3 / 4 n/y/n), it can be easily shown 
that A and A 1 is 0(r). 

Now, nodes that are in S (i.e. physically inside S) can only be connected to those nodes in S c that lie in A 1 . 
Similarly, nodes in S° that are connected to nodes S must lie in A . Thus, nodes that can communicate across 
the cut (S,S C ) must he within a region of area 0(r). For the protocol model, if a node transmits, nodes within 
distance r(l + 77) of the receiver must not transmit. That is, each transmission effectively silences nodes within an 
area of ©(r 2 ). Thus, at any given time, the maximum number of simultaneous transmissions between S and 5 C is 
G(l/r). This along with \S\, \S C \ = 9(n) implies that 

* < qt(s) = ()[ — ) O 



0(re 2 ) \n 2 rj log 3 / 4 n, 

For the lower bound, we identify a grid subgraph of G{n,r) with r = 0(log 3 / 4 n/^/n). Consider a grid graph 
G n of yjn x y/n nodes with each node connected to one of its four neighbors (with suitable modifications at the 
boundaries). The nodes of G n are placed in a uniform manner in a unit square; each node is at a distance \ j\fn 
from its neighbors. Now consider a minimax matching between nodes of G n and re randomly placed nodes in the 
unit square, where a minimax matching is a perfect matching between the n nodes of G n and the nodes of G(n, r) 
with maximum length minimized. Leighton and Shor [21] established that the maximum edge length in a minimax 
matching, say r*, is G (log 3 / 4 nj^fn) with probability at least 1 — 1/n 4 . Now we identify the subgraph G' n (with 
grid graph structure) of G(n, r) as follows. G' n has all n nodes. Consider the minimax matching between G n and 
G(n,r). If a node of G(n,r) is connected to node number m of G n , then renumber it as rre to obtain nodes of 
G' n . Now by setting r > r* + clearly a node m and m' are connected in G' n if they are connected in G n . 

Thus, we have established that G n C G' n . Now, we will focus only on the edges of G(n, r) that belong to G' n 
and provide them with positive capacity by an appropriate communication scheme that is feasible for the restricted 
protocol model. For this, note that in G(n, r) each node is connected to at most 0(log 3 / 2 n) nodes with probability 
at least 1 — 1/n 4 (using Chernoff's bound and Union bound) for large enough n. Hence, using a simple TDMA 
scheme based on vertex coloring of G(n,r), each node gets to transmit once in every 6(1/ log 3 ' 2 n) time slots. 
This transmission can be along any outgoing edge. Since, we are interested in providing positive capacity to only 
at most 4 outgoing edges, we have established that there is a simple TDMA scheme which provides 0(1/ log 3 / 2 n) 
capacity to each edge of a grid subgraph of G(n,r). To complete the proof, we recall that the conductance of a 
grid graph is Q(l/^/n) [22]. That is, 

^tn \ ■ ^ieSjeS" 1 {(i,j)eE} 
<£>(G n ) = mm , = 



s \S\\S C \ " V^ 3/2 

Now, putting all the above discussion together we have the following. 



VP = sup min 

c&cscv 



cecscv \S\\S- 
1 



log 3 / 2 n 



(re 3 / 2 log 3 / 2 n) ^ ^ 

In summary, upper and lower bound on ^ along with Theorem [3] implies the Lemma [6l ■ 

Now, we discuss briefly delay. In [4], delay was defined as the average number of hops per packet, and the packet 
size was assumed to scale to an arbitrarily small value. For any communication scheme feasible for the protocol 
model with maximum transmission radius r = ©(log 3 / 4 n/^/n), the maximum number of transmissions per unit 
time is upper bounded as 0(n/log 3 / 2 n). Using this and Theorem |4]we obtain the following result immediately. 
Lemma 7: The delay D{n) for any scheme achieving f* = Q ^ n3/2 j~p72 n ) is bounded above as 

D{n) = O (\/nlogn) . 



13 



IV. Gaussian Fading Channel Model 

In the previous section, we assumed that the wireless network was defined by two graphs G and G D . We 
extended the results of Leighton and Rao to wireless networks modeled by a combinatorial interference model; 
this mainly exploited the fact that all possible transmission schemes could be described in terms of routing over a 
set of capacitated graphs, where the set of edge capacity vectors belonged to the convex hull of a finite number of 
vectors. Thus, in this sense, the inherent discrete nature of the model worked to our advantage. 

While the combinatorial interference model can allow for arbitrary scheduling and routing schemes, it does not 
model all the degrees of freedom in a wireless network. Specifically, the results are not information theoretic. In 
this section, we provide an information theoretic characterization of PMF in a wireless network with Gaussian 
fading channels. The techniques for the combinatorial model can be easily extended to obtain a feasible scheme 
and a lower bound on the maximum PMF /*. However, for information theoretic upper bounds we have to work 
harder, especially to obtain a bound that relates to the lower bound and allows us to quantify the gap. 

Our key contribution is in quantifying the suboptimality of the UMF/PMF for a simple feasible scheme and 
an upper bound on the UMF/PMF for an arbitrary network topology, in terms of a simple graph property. The 
bound is general when channel side information (CSI) is assumed to be available only at the receiver. For AWGN 
channels, we quantify only for UMF, and when the SNR is low enough. To the best of our knowledge, this is the 
first such result which guarantees that a feasible scheme achieves rates within a certain factor of an outer bound for 
an arbitrary graph. We also illustrate these results through applications. The results hence obtained are interesting 
in their own right. 

Our main approach is as follows. We construct two directed capacitated graphs G u and G L for the given wireless 
network. The graph G u is such that the capacity (defined appropriately later) of each cut in G u upper bounds the 
corresponding cut-capacity in the wireless network. The graph G L is such that there exists a communication scheme 
that simultaneously achieves the capacity of each edge in G L , and the ratio of capacity of each cut in G u and 
G L is bounded above by a quantifiable term. This leads to an approximate characterization of PMF in an arbitrary 
wireless network with Gaussian fading channels. Moreover, the feasible scheme that induces the capacities in G L 
supports PMF which is within a quantifiable factor of the optimal. 

A. Channel Model 

This is similar to the model in, for example, [23]. We have V = {1,... ,n} wireless nodes with transceiver 
capabilities located arbitrarily in a plane. Node transmissions happen at discrete times, t E Z+. Let Xj(i) be 
the signal transmitted by node i at time t G Z + . We assume that each node has a power constrain^ such that 
limsup^v^oo i Ylt=l \Xf(t)\ < P. Then Yi(t), the signal received by node i at time t, is given by 

Y i (t) = Y lk ^H ik X k {t) + Z i (t), (5) 

where Zi(t) denotes a complex zero mean white Gaussian noise process with independent real and imaginary parts 
with variance 1/2 such that Z^it) are i.i.d. across all i. Let denote the distance between nodes i and j. Let 
Hik(t) be such that 

H ik {t) = \Zg(r ik )H ik (t), 

where Hi k (t) is a stationary and ergodic zero mean complex Gaussian process with independent real and imaginary 
parts (with variance 1/2). It models channel fluctuations due to frequency flat fading. Also, g(-) is a monotonically 
decreasing function that models path loss with g(x) < 1 for all x > 0. We assume also that the iJjfc(t)'s are 
independent. 

B. Graph Definitions 

Consider the following two graphs induced by a wireless network of n nodes: 
(1) K n is the fully connected graph with node set V; 

7 For notational simplicity we assume that each node has the same power constraint. The general case, where each node has different 
maximum average power can be handled using identical techniques. 
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(2) G r is the graph where each node i G V is connected to all nodes that are within a distance r of i. Let E, 
denote the edge set of G r . Let A(r) be the maximum vertex degree of G r . Finally, define 

r* = min{r : G r is connected}. 



C. Preliminaries 

In the analysis in this section, we utilize the following two simple lemmas. 
Lemma 8: Given Xj G (0, 1), 1 <i<N, 

Eli log(l + V^) < V^N^Eli^gil + x,). 
Proof: For any x G (0, 1), x/2 < log(l + x) < x, so 

EiIiiog(i + ^) < EliV^< ^vTeEi^ (6) 

where ([6]) follows from Cauchy-Schwarz inequality. ■ 

Lemma 9: For any x > 0, a G (0, 1), ^ log(l + ax) > log(l + x). 

Proof: Define f(x) = ± log(l + «x) - log(l + x). Note that /'(x) > for x > and f(0) = 0. ■ 

D. Results 

We obtain bounds on the maximum PMF for three different cases: 

(1) fading channel with AWGN, and channel side information (CSI) available only at the receiver, 

(2) deterministic (no fading) additive white Gaussian noise (AWGN) channel, and 

(3) fading channel with AWGN, and CSI available at both the transmitter and the receiver. 

The exact bounds for the above cases are different, but the analysis and bounding techniques are similar. 



1 ) Random Fading with Rx-only CSI: We first obtain bounds on the PMF for Gaussian channels with random 
fading under the assumption that CSI is available at the receiver, but not the transmitter. We then relate the bounds 
for PMF, and show that the gap can be quantified well, and under very general assumptions. We note that this is 
the case for which we can obtain the strongest results. 

Theorem 5: With channel state information (CSI) only at receivers, /* is bounded as follows: 

£^ e s^(iog(i + J p|^| 2 )) 



ft < min 

Jn ~ scv 



tt(S)it(S c ) 



sup 

r>r*, r/>0 



1 + A(r)A(r(l +r])) 



EieSjeS»%fle^ Elo g( 1 + T 



mm 

scv 



mill 



-nPg(r(l+n)) 



logp 7r 7r(5)7r(5 c 



Theorem [5] provides bounds on /* which relates to the "cut capacity" of appropriate capacitated graphs. Specifi- 
cally, we can compute the information theoretic upper bound (for any PMF) in polynomial time using flow arguments, 
and by solving an LP as detailed in Sec. \IV-E\ However, it is not clear how tight these bounds are. We now quantify 
the gap between the upper and lower bounds. 

Corollary 2: For any r > r* , denote S(r) = maxj Ylj-r >r ^9( r ij)- Then, 

"I- T ^ ^ </;<(l + 7(0)T, 



(l + A 2 (r))(l + <5(r))log^ 



where 



T 



mm 

scv 



E^ e ^:,, 3 <^lQg(l + P|^| 2 ) 

ir(S)ir(S c ) 
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and 

IWs.:r y >rElog(l + P|flii| 2 ) 

7(rj = max 



5cy:7r(5) j7 r(5C)>0 £ ieSJeS c :r . .< r E log(l + P|H;j | 2 ) ' 

Note that both 7(r) and d(r) are decreasing functions of r, while A(r) is an increasing function of r. Also, 
since power typically decays as l/r a for 2 < a < 6, while for uniformly distributed networks A(r) grows only 
linearly with r, the decay of j(r) and S(r) is much faster than the growth of A(r). Hence, for r large enough the 
gap is dominated by the term logp^l + A(r) 2 ). Specifically, assume that there exists an e > such that the graph 
G e = (V, E e ) is connected, where E e = : E [log (l + P\H i:j \ 2 )] > n~ e / 2 }. Then the above bound for UMF 

reduces to [24] 

H^)i^j =/ =0(min - cutR) ' 

E ieS ,eS-r <r El0g(l + P|H i3 | 2 ) 1 

where min-CUiR = minscv : ^tw^ct , and r € is such that S(r e ) < ^t+t- 

Proof: [Theorem [5J We first prove the upper bound. Following the steps in the proof of Theorem 2. 1 in [23] 
and using (1 + Yl7=l a «) — EKLi(l + a i) f° r a i > 0> we obtain that for A £ A, 

EieSjeS^ij < n max E[\og det(I + H S Q S H* S )] 

< E, eW E(log(l + P|^| 2 )). (7) 

Now, for any PMF M = M(/ 7r ,vr), it must be that Eie5,je5 c M u = fn7r(S)ir(S c ). Hence, for any such PMF 
M(f, 7r) G A, the upper bound in the Theorem holds. 

To establish the lower bound, we construct a transmission scheme for which the PMF is greater than or equal 
to that in the lower bound. For r > r* , consider the graph G r = (V,E r ) on the n nodes defined above. We use 
A(r(l +77)) to denote the maximum vertex degree of the graph G r (i +ri y Now, consider the following transmission 
scheme. A node % can transmit to a node j only if rij < r. Also, when a node i transmits, no node within a 
distance r(l + 77) of the receiver can transmit. Thus, when a link G E r is active, at most A(r(l + 77)) nodes 
are constrained to remain silent, i.e., at most a = A(r(l + r/))A(r) links are constrained to remain inactive. Hence, 
the chromatic number of the dual graph is at most (1 + A(r(l + rj))A(r). In addition, we assume that the signal 
transmitted by each node has a Gaussian distribution. For any given link that transmits data at a particular time, 
we treat all other simultaneous transmissions in the network as interference. Now focus on any one link, say link 
(1,2) between node 1 and 2, without loss of generality. We claim the following. 

Lemma 10: For the above scheme, the following rate on link (1,2) is achievable: 

A 12 = «- 1 Elog(l + IT ^g^ ). 
We prove Lemma [10] later. First we explain how it implies the proof of Theorem |5J A similar analysis holds for 
other links that (1, 2) in E r . Thus, for graph G r the following rate are achievable on link G E r ■: 



a 



1 ( P\Hd 2 \ 

^Elog 1 + 7T 

6 V l + nPg(r(l+r]))J 



Now given the capacitated graph G r , we can use classical wireline network based routing algorithms for obtaining 
a product multicommodity flow that is lower bounded by the following quantity: 



fLB(r,v) = ^ 



1 



1 + A(r)A(r(l + 77)) 
This implies the following lower bound on /* 



mm 



T,ies,jeS° %j)e£,. Elo g ( 1 + l +n P^r(l+r,)) 



scv logp 7T TT(S)ir(S c 



K> sup /lb (r, 77). 

r>r* ,r]>0 

This is precisely the claimed lower bound in the statement of Theorem [5j and thus completing the proof. 
Proof: [Lemma [TOl We will use the following result, that follows directly from Theorem 1 in [25]. 
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Theorem 6: Consider a complex scalar channel where the output Y when X is transmitted is given by 

Y = hX + Z + S, 

where Z is a complex circularly symmetric Gaussian random variable with unit variance, and S satisfies E[5*5] < 

P. Also, h is zero mean and i.i.d over channel uses. If X is a complex zero mean circularly symmetric Gaussian 

P\h\ 2 ' 



random variable with E[X*X] = P, then I(X; (Y, h)) > Elog (1 + j 

We consider a transmission scheme where the signal transmitted over each link, when active, is a complex zero 
mean white circularly symmetric Gaussian with variance P. Moreover, we assume that the transmissions on all 
links are mutually independent. Let ti,t 2 , ■ ■ ■ denote times at which link (1,2) is scheduled. Hence, at any such 
time t G {ti, t 2 , ■ ■ .}, the received signal at node 2 is given by 

Y 2 (t) = H 21 (t)X 1 {t) + E W H 2k(t)X k (t) + Z 2 (t). 

Using the mutual independence of transmissions and zero mean property along with the construction of the 
scheduling scheme, 

E\E W H 2k (t)X k (t) + Z 2 (t)\ 2 <l+nPg(r(l+ri)). 

From Theorem [6l 

J(Xi(t); (Y 2 (t),H 21 (t))) > Elog (l + {1+n P ^ +v)) ) . (8) 

Since the channel is assumed to be i.i.d. over channel uses, a random coding argument can be used to achieve this 
rate with a probability of error that goes to zero as the block length goes to infinity. 

Combining this with the time-sharing between different sets of links described above, since each link gets to 
transmit at least once in a times slots, or at least 1/a fraction of the time, it follows that 

/ P\H I 2 

A12 > a _1 Elog 1 + ! 2J| 



1 + nPg(r(l + r?)) 



Proof: [Corollary 13 Consider any S such that tt(S),tt(S c ) > 0. Then, 

Cut(S, S c ) = EieSjeS* E (log(l + P\H^)) . 
< (l+7(0)E, eW :n J <rE(Ml + J P|^| 2 )) , (9) 

where the second line follows from the concavity of the log function, Jensen's inequality, log(l + x) < x for x > 
and definition of 7(7"). Thus, 

The upper bound then follows from the upper bound in Theorem [5] 

Next, we consider the transmission scheme that led to the lower bound in ([8]) with 77 = 0. Note that in ([8]>, 
we used the term nPg(r(l + 77)) as a bound on the interference power. However, here we consider the actual 
interference = X^fcev-r k >r P9i r jk) f° r a transmission from i to j. Note that Ijj < 6(r). Now, by Lemma |9j 
we have 

E [log (l + ^l^)] > t^IE (log(l + .P|^| 2 )) • (ID 

Using (© and (TTTb along with the lower bound obtained via time-division scheme that led to ([8]), the lower bound 
in Theorem [5] gives us 

j = il mm * 



scv tt(S)it(S c )(1 + A(r) 2 )(l + 5(r)) logp^ 
T 

(1 + A(r) 2 )(l +5(r))logp„ 



, -, N9W T , Tl — ). (12) 
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2) Deterministic AWGN Channels: We now consider an AWGN channel without fading, i.e., we have H^j = 1 
w.p. 1, Vfc, j = 1, ... ,n. We first obtain the following set of bounds on maximum PMF using standard arguments. 
Theorem 7: The maximum PMF /* is bounded as follows. 



sup 

r>r* , r/>0 



f* < min 
j7T - scv 



1 



2E 



l0g(l + y/Pjfcj)) 



1 + A(r)A(r(l + t?)) 



ir(S)Tr(S c ) 



mm 



1 _|_ p g( r y) 

1 l+nP S (r(l+r/)) 



logp n ir(S)7r(S c 



Next, we present a Corollary of Theorem [7] which characterizes the tightness of the above bound for UMF for 
low signal to noise ratio (SNR). 

Corollary 3: Define J(r) = min{7 > : Ylj-r >r P9i r ij) — ^> f° r au r (^) = min{r > : J(r) < 5} for 
5 > 0. Then, 



n 



n 



(1 + <y)A(r(5))(l + A(r(<5)) 2 ) logn 



< /* < 2T + O - 
n 



where 



T 



mm 

Ucv 



£ietfjetf«:r y <r(a) ^SC 1 + V /P 5(^j)) 



IW C I 

We now present the proofs of Theorem |7] and Corollary [3] The main idea in the proof of Theorem [7] is to neglect 
interference to upper bound achievable rates on links, and to construct a transmission scheme to induce achievable 
rates on the links. In particular the scheme that we construct consists of time sharing between multiple transmission 
schemes, each of which enables direct transmissions between nodes that are separated by at most distance r. Then 
the lower bound on /* is obtained by routing over graph G r , where each edge has a capacity given by this time 
division scheme. 

Proof: [Theorem 13 We first prove the upper bound. In order to bound the sum-rate across each given cut, we 
refer to the proof of the max-flow min-cut lemma in [8], which yields for any S C V and A £ A, 



£ J6 



where Xj = E ieS VdFj 



Xj . We therefore deduce that 



Aii<£ i6S . log(l + E(|X,f)), 



E 



< Ejes- log[l + EiMS V»M |EpQX fc 

< E, 6 s c iog[i + P(E ? : e sVsM) 2 ], 



since \E(XiX k )\ < ^/P~P\ < P. Finally, we obtain 

Eiesjes* X ij ^ EieS°21og(l + V^Eies V^i)) 

Now, for any PMF M = M(/ 7r ,vr), it must be that Ei 6 s,j£S° M y = f^iS)^ )- Hence, for any such PMF 
M(f, it) G A, the upper bound in the Theorem holds. 

To establish the lower bound, we construct a transmission scheme for which the PMF is greater than or equal 
to that in the lower bound. For r > r* , consider the graph G r = (V,E r ) on the n nodes defined above. We use 
A(r(l + 77)) to denote the maximum vertex degree of the graph G r (i +rj y Now, consider the following transmission 
scheme. A node i can transmit to a node j only if rij < r. Also, when a node i transmits, no node within a 
distance r(l + 7/) of the receiver can transmit. Thus, when a link (i, j) G E r is active, at most A(r(l + if)) nodes 
are constrained to remain silent, i.e., at most A(r(l + 7/))A(r) links are constrained to remain inactive. Hence, 
the chromatic number of the dual graph is at most (1 + A(r(l + 7?))A(r). In addition, we assume that the signal 
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transmitted by each node has a Gaussian distribution. Then, subject to the maximum average power constraint, for 
any node pair i,j, such that < r, the following rate is achievable from i — > j: 

l+nPg(r'(l+^))) 

l + A(r)A(r(l + T/)) " U } 

Note that the interference is due to at most n nodes and all the interfering nodes are at least a distance r(l + 77) 
away from the receiver. We now consider routing over the graph G r , where each edge has capacity Ay. The 
lower bound then follows from the lower bound in Theorem |2] ■ 



Proof: [Corollary O Consider any cut defined by (S, S c ). Due to the symmetry of the upper bound in Theorem|7J 
without loss of generality, assume \S\ < n/2. Consider any 5 such that r(6) > r*. Then, 



Cut(5, S c ) = Eiesjes* Ml + s/PgJ^)) 
= Eies,jes-.r u <r(S) lo § (! + V p 9( r ij)) + E l esjes^>r{S) lo § i 1 + V p 9( r ij)) 

log (1 + y/P^J) + \S\S, 



(14) 



where the last step follows from the definition of r(S), and log(l + y/x) < x for x < 1. Hence, the upper bound 
in the Corollary follows from the upper bound in Theorem [7] Since we assume Pg{rij) < 1 for all i and j, from 
Lemma [U we have 



ieS,jeS":r ij <r(S) 



1 + \Pg(rij) 



< 



l2A(r(5))\S\ log(l + P 5 (^)). 

ieS.,]eS":r z:1 <r(S) 



(15) 



For the lower bound, consider the choice of r = r(<5) and r\ = for the scheme described in the proof of 
Theorem |7J Then, the interference during data transmission from i to j, Iy = YlkeV-r k >r(8) Pd( r jk) — Now, 
Lemma |9j implies that 



log 1 + 



Pg(r. 



1.1. 



> 



1 



l + Iij J ~ 1 + 5 



log(l + Pg( rij )). 



(16) 



Using an appropriately modified lower bound in Theorem [7] for the choice of r = r{5), r\ = 0, it follows that 

£ie£WeSe:r 4 ,<r(i) l0g(1 + P 9( r ^)) 



r = 



sgv (l + 5)(l + A(r((J)) 2 )logn|5||5 c 
T 2 n 



(l + ^)A(r(5))(l + A(r(«5)) 2 )lognJ ' 

where the second step follows from (fT5l) . The lower bound in Theorem |7] then implies the lower bound in the 
Corollary. This completes the proof. ■ 



3) Random Fading with CSI at both Tx and Rx: We now obtain bounds on the PMF for a Gaussian channel 
with random fading when CSI is available at both the transmitter and the receiver. Qualitatively, these bounds are 
very similar to the case of deterministic AWGN channels. The main result is as follows. 

Theorem 8: With CSI at both transmitters and receivers, /* is bounded as follows. 

. Eiesje3^0osO- + VP\H ji \) 

t < mm . 

J * ~ scv ir(S)ir(S c ) 

The lower bound for the receiver only CSI case is a (weak) lower bound for this case as well. 



Proof: The upper bound follows again from the proof of Theorem 2.1 in [23], from which we deduce that 
for any A £ A, 

EieS,jeS<= A *i ^ E [ max Qbo,Q„<P log det(J + H s QsH* s )] 
< J2jeS' E [ max QbO,Q l ,<plog(l + hjQsh*)}, 
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where hj is the j th row of H. Since hjQsh* is maximum when (Qs)ik = P f° r ai l i,k E S, we obtain, following 
the steps of the proof of Theorem |7J 



>:, a a,, < E, e 5^io g (i+p(E, e5 i^i) 2 ] 

< E^ 2E1 °g(i + V^E i6S l^il)] 

< Ei 6S , ie se2E(log(l + VP|^|), 



so the upper bound on /* follows from Theorem [2] ■ 



E. Computational Methods 

We discuss the implications of the bounds for the case of CSI availability at the receiver only as stated in 
Corollary |2] Similar implications follow for the case where CSI is available to both transmitters and receivers as 
well. 

Corollary |2] shows that an upper bound on /* can be obtained via the maximum PMF on graph K n , where each 
edge has a capacity log(l + P\Hij\ 2 ), and there is no interference; specifically, logn times the PMF thus 
computed for K n is an upper bound on /*. The lower bound is obtained via routing on G r with edge having 
capacity n^MH^^+s) • Hence, the PMF on G>(<j) is a lower bound on /*. Both the above computations can be 
done by solving an LP in polynomial time. Moreover, the ratio of the bounds is quantified in |2] We note that 
such an approximation ratio could be obtained easily for the combinatorial interference model using node coloring 
arguments. The arguments here are more complicated, as detailed in the proof of Corollary |2] 

F. Application 

We now apply the information theoretic characterization of PMF in the previous subsection to obtain a scaling 
law for average UMF in a geometric random network with a fading channel, and when CSI is available at the 
receivers. The scaling law we obtain is along similar lines to those that exist in the literature. Similar bounds can be 
obtained when CSI is available both at the transmitter and the receiver or when the channels are AWGN channels. 

We consider a geometric random graph model with a constant node density: n nodes are placed uniformly at 
random in a torus of area n (and not unit area). Thus the distance between two nodes is a random variable taking 
values in (0, 0(y / ra)). We assume that all nodes have the same transmission power equal to 1, i.e. Pj = 1 for all 
1 < i < n. We state the following result characterizing /*. 

Lemma 11: Consider the Gaussian channel model with random fading and CSI available only at the receivers. 
Let g(r) = (1 + r)~ a , a > 3 and P = 1. Then for a geometric random graph with constant node density (described 
above), the average (over random position of nodes) /* is bounded as 



if Pr(\Hij\ 2 > 0) > 7 for some strictily positive constants j3, 7 (independent of n) for all 1 < i,j < n. (Note that 
the condition is for the normalized channel gains HijS and not the actual gains HijS.) 

Proof: We use Theorem [5] to evaluate the bounds. First we obtain an upper bound by evaluating the bound 
of Theorem [5] for a specific cut {U, U c ). Then, we evaluate lower bound by relating it to an appropriate grid-graph 
as in Lemma [6] 

Now, we consider the upper bound. Consider a horizontal line dividing the square of area n into equal halves. 
Let U be set of nodes that lie in bottom half, and so U c is the set of nodes that lie in the top half. From Theorem 
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|5J we have 

r\u\\u c \ < Yl emi+pi^ 12 ' 



< Yl io g (i+™[i^i 2 ]) 
= J2 io g (i+p 5 (^)) 

ieu,jeu° 

= Yl^+^y^am 1 ^}, (is) 

where we have used Jensen's inequality, log(l + x) < x for all x > 0, and the hypothesis of Lemma. Since, the 
nodes are thrown uniformly at random, the expectation of each term in ( fT8l ) for a pair is the same. Using 
linearity of expectation, we obtain that 

EE( 1 + %)~ ai {iet/}l{i6^}] = n(n-l)E[(l + r 12 )- a l { i el / } l {2el /c } ] 



a sds y/ndr 




= O(^), (19) 

where we used the fact that for a > 3 the last integral is bounded above by a constant. The above evaluation can 
be justified as follows. First note that Pr(l G U, 2 £ U c ) = 1/4. Given { 1 G U,2 G U c }, node 1 in the bottom 
rectangle and node 2 in the top rectangle are uniformly distributed. Now, consider a thin horizontal strip of width 
dr and length y/n at distance r below the horizontal line dividing the square (and inducing U, U c ). The node 1 G U 
belongs to this strip with probability 2y/ndr/n. Now, the node 2 is at distance at least r from node 1. Consider a 
ring of width ds, centered at node l's position and of radius s > r. The area of this ring is 2irsds. The probability 
of node 2 being in this ring is bounded above by Ansds/n. When the above described condition is true, the nodes 
1 and 2 are at distance s. Integrating over the appropriate ranges justifies the final outcome ( fT9l ). 

Now, it is easy to see that under any configuration of nodes, /* = 0(n 2 ) since g(r) < 1 for any r > 0, P = 1 
and elementary arguments. Let event A = {\U\\U C \ = 0(n 2 )}. Using Chernoff's bound, it is easy to see that (with 
appropriate selection of constants in definition of A) for large enough n, we have 

Pr(A) = 1 - 1/n 6 . 

Using this estimate and bound /* = 0{n 2 ) we obtain that 

E[/*] = E[f*l A ]+E[tl A o] 

< E[ri A ]+o(- A 



n 2 ) In 4 



_ e| -E[/win + Q (l 



n 2 J \n 4 
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Next, we prove the lower bound. For this we construct a graph with achievable link capacities for which the 
average /* is lower bounded as claimed in the Lemma. Consider r = B(logn). Then the corresponding G r , which 
is the geometric random graph G(n, r), is connected with high probability (at least 1 — 1/n 4 by appropriate choice 
of constants in selection of r). For this choice of r, using the Chernoff and Union bounds it follows that with 
probability at least 1 — 1/n 4 , 

A(2r) = 6 (log 2 n) . 

Again, we can identify a grid graph structure as a subgraph structure of G r based on the argument used in Lemma 
[6l Denote the edges of this grid sub-graph structure as E. We note that 0(1) edges are incident on each of the n 
nodes that belong to E (which is a property of the grid-graph structure). Next, we design a feasible transmission 
scheme for which each edge in E can support a transmission rate of f2(log _a n). 

Specifically, we consider a TDMA schedule for the graph G r similar to that described in the proof of the lower 
bound for Theorem[5] It is easy to see that G2r can be vertex colored using 0(A(2r)) colors. We use a randomized 
scheme to do TDMA scheduling as follows: in each time-slot, each node becomes tentatively active with probability 
1/A(2r) and remains inactive otherwise. If a node becomes tentatively active and none of its neighbors in G2 T is 
tentatively-active, then it will become active. Else, it becomes inactive. All active nodes transmit in the time-slot 
simultaneously. It is easy to see that each node transmits for 0(1/A(2r)) fraction of the time on an average. The 
randomization here is used to facilitate the computation of a simple bound on the average interference experienced 
by a node due to transmissions by nodes that are not its neighbor. 

Now under the above vertex coloring, each node gets to transmit once in 0(A(2r)) time-slots on average at 
power 0(PA(2r)) = 0(A(2r)). We wish to concentrate on transmissions for edges that belong to E, which are 
a subset of edges of G r . For any such transmission, say from u — > v with (u, v) £ E, according to the above 
coloring of G2 r no other node within distance r of v is transmits simultaneously. Also, any node that is at least 
distance r away from v can be active with probability at most 1/A(2r). Hence, the average power corresponding 
to the interference received by node v, say I v , can be bounded above as follows: 

\j:r vj >r ) \j:r VJ >r ) 

where we used the fact that each node transmits at power PA(2r) = A(2r) for 1/A(2r) fraction of the time and 
E[|iZjj| 2 ] = g{?ij). By another application of Chernoff's bound and union bound, it can be shown that the number 
of nodes in an annulus around node v with unit width and radius R for R G N, r < R < y/n, is Q(R) with 
probability at least 1 — 1/n 5 . Then it follows that 




S(r*j) = V g(R)R] = O (,--»+ 2 ) , 



£ 

j:r vj >r 

where we have used fact that a > 3. Using this in (12T1 ). we obtain that with probability at least 1 — 1/n 5 / 2 we have 

I v = O ( r - a+2 ) = O (log 2 - a n) . (22) 

That is, I v — ► as n — ► oo for a > 3. Thus, by selection of large enough n, I v can be made as small as possible. 
That is, when transmission from u — > v happens, the average noise received by node v due to other simultaneous 
transmission is very small, say less than 5 for some small enough 5 > 0. 

Given this, the arguments used in Lemma \10\ imply that when u transmits to v at power PA(2r) once in 
0(A(2r)) time-slot, considering other transmissions as noise, we obtain that the effective rate between u — > v is 
lower bounded as 

PA{2r)g(r uv )\H uv \ 2 



X u ^ v = Q A(2r)" 1 E 



log 1 + 



l + Iv 



Now, H uv is independent of everything else and 

Pr(|^™| 2 > P) > 7, 
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for some positive constants f3,j as per our hypothesis. Therefore, use of Lemma [9] implies 

A_>^fA(2^1ogfl + PA(2r) ^ 



l + Iv 

Further, I v < 5 for small enough 5. Therefore, another use of Lemma [9] implies that 

X u ^ v > n (A(2r)- X log (1 + PA(2r)g(r uv ))) . 

Now 

Pg(r uv )A(2r) = ft (log 2 ~ a n) , 

where we have used the fact that r uv < r and g(-) is monotonically decreasing. For x £ (0, 0.5), log(l + x) > x/2. 
Therefore, for a > 3 

A u ^ = ft(A(2r)- 1 log 2 - a n). 

Since A(2r) = G(log 2 n), we have established that the effective capacity of transmissions for each edge under the 
above described TDMA scheme is ft(log _Q n). That is, each edge of E gets capacity at least Q(log~ a n). Now 
recall that a grid graph with unit capacity has /* lower bounded as C^T^togrl ) • Hence, using this routing of 
UMF along edges of E with capacity f2(log~ a n) we obtain that 

r = n- 1 



n 3 / 2 log 1+a n, 

By careful accounting of probability of relevant events above and Union bound of events will imply that the above 
stated lower bound on /* holds with probability at least 1 — 1/n 2 . Since /* > with probability 1, it immediately 
implies the desired lower bound of Lemma: 

1 



e[/*] = n 

This completes the proof of Lemma [TTJ 



n 3 / 2 log 1+Q n 
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